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We show in detail how Richardson's exact solution of a discrete-state BCS (DBCS) model can be 
recovered as a special case of an algebraic Bethe Ansatz solution of the inhomogeneous XXX vertex 
model with twisted boundary conditions: by implementing the twist using Sklyanin's K-matrix 
construction and taking the quasiclassical limit, one obtains a complete set of conserved quantities, 
Hi, from which the DBCS Hamiltonian can be constructed as a second order polynomial. The 
eigenvalues and eigenstates of the Hi (which reduce to the Gaudin Hamiltonians in the limit of 
infinitely strong coupling) are exactly known in terms of a set of parameters determined by a set 
of on-shell Bethe Ansatz equations, which reproduce Richardson's equations for these parameters. 
We thus clarify that the integrability of the DBCS model is a special case of the integrability of the 
twisted inhomogeneous XXX vertex model. Furthermore, by considering the twisted inhomogeneous 
XXZ model and/or choosing a generic polynomial of the HiS as Hamiltonian, more general exactly 
solvable models can be constructed. - To make the paper accessible to readers that are not Bethe 
Ansatz experts, the introductory sections include a self-contained review of those of its feature which 
are needed here. 



I. INTRODUCTION AND SUMMARY 

In a series of pioneering experiments in the midfjl990's, Ralph, Black and Tinkhani observed a spectroscopic gap 
indicative of pairing correlations in Al nanograinsB that were so small that their electronic excitation spectra were 
discrete. Their results inspired a growing number of theoretical studies of superconducting pairing correlations in 
nanograins with a fixed number electrons (see Refs. |^,^ for recent reviews). These works are based on a model, to 
be called discrete-state BCS model (DBCS model) below, described by a reduced BCS Hamiltonian for a discrete 
set of doubly-degenerate energy levels, with a pairing interaction that scatters pairs of electrons from one level to 
the next. The DBCS model was solved exactly by Richardson in a series of papers starting in 1963u he explicitly 
constructed all eigenstates and eigenenergies of the DBCS Hamiltonian in terms of a set of energy parameters whose 
values are determined by (numerically) solving a set of algebraic equations, to be called "Richardson's equations". 
Though his work had, for a long time, been overlooked by the condensed matter community, it has recently received 
increasing attention in the context of studying pairing correlations in nanoscale Al grains, where the existence of an 
exact solution has turned out to be as useful as it had been unexpected. 

The existence of an exact solution to a nontrivial model of course immediately raises the question whether it is 
related to any of the standard ways of exactly solving solvable models. The goal of this paper is to show that this is 
indeed the case: Richardson's solution of the DBCS model is a special case of an algebraic Bethe- Ansatz solution of 
the so-called inhomogeneous XXX vertex model with twisted boundary conditions. 

This insight builds upon a seoes of recent observations regarding exact properties of the DBCS model: In 1997, 
^ . Cambiaggio, Rivas and SaraccnocI showed (though unaware of Richardson's work) that the DBCS model was integrable, 
■ and explicitly constructed all the constants of the motion [cf. Eq. ( p^ ) below]. In 2000, Amico, Falci and FazioQ realized 
rN ' that the DBCS integrals of motion are in fact very similar to the integrals of motion of the XXX Gaudin model [cf . ( |Tl| ) 
below], differing from the latter only by an additional Sz term, and that Richardson's equations are very similar to 
■ ■ ■ ' the so-called Gaudin equations, differing from the latter only by an additional constant term. Now, it has long been 
known (see, e.g., chapter 13.2 of Gaudin's bookCI) that the Gaudin model can be derived from the inhomogeneous 
XXX vertex (IXXX) model with periodic boundary conditions, by taking the so-called quasiclassical limit, and that. 
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correspondingly, the Gaudin equations can be derived by taking the quasiclassical limit of Bethe Ansatz equations of 
the IXXX model. Since-iRichardson's Ansatz satisfies the Gaudin equations modified by the additional constant term, 
Amico, Falci and Faziocl referred to Richardson's Ansatz as an "off-shell Bethe Ansatz" , i.e. an Ansatz not satisfying 
the Bethe-equations of the original XXX model, but of a modified version thereof. (The off-shell Bethe Ansatz was 
originally introduced by Babujian and Flume in a context quite different than finding cigcnstatcs and eigenvalues of 
integrai models, namely, to solve Knizhnik-Zamolodchikov differential equations arising in conformal Wess-Zumino 
modelfS.) 

In this paper, we address the following question: can one construct a vertex model, integrable by the algebraic Bethe 
Ansatz (ABA), whose quasiclassical limit directly gives the DBCS model, in other words, which is directly solved by 
a normal "on-shell" Bethe Ansatz? The answer is positive: we show that the sought-after model is an IXXX model 
with twisted (instead of periodic) boundary conditions, which we shall call the TIXXX model; its transfer matrix 
yields, in the quasiclassical limit, a complete set of conserved quantities. Hi, from which the DBCS Hamiltonian can 
be constructed as a second order polynomial. Our emphasis on twisted boundary conditions is the main difference 
between our work and that of references |^J^. We implement the twist using the boundaj:K-.K-matrix construction 
of Sklyanin, which he introduced while developing his method of separation of variables, L2rli3 an alternative (and in 
some cases more powerful) way to the ABA for constructing wave functions. In fact, Sklyanin himself mentioned in 
a side remark in Ref. |ll| that the quasiclassical limit of the IXXX model with twisted boundary conditions (using a 
diagonal iiT-matrix) produces a modified version of the Gaudin model (though he was not aware, at the time, of the 
connection of the latter to the DBCS model). 

We hope that our work fully clarifies the origin of the integrability of the DBCS model by explicitly constructing the 
integrable TIXXX model from which the latter can be derived. Moreover, by this construction we pave the way for 
using the powerful algebraic Bethe Ansatz machinery to calculate various quantities that have not yet been studied for 
the DBCS model. For example, there has recently been great progress in using the ABA to calculate matrix elements 
(or form factors) and correlation fupctions in vertex models, which, by building upon our work, could now fruitfully 
be apphed to the DBCS model, tooll3. 

Our work also suggests ways for constructing integrable generalizations of the DBCS model, by considering other 
vertex models with twisted boundary conditions. In fact, one such gencpjalization, which is Bethe- Ansatz solvable, has 
recently been found independently by Amico, Di Lorenzo and Osterlotcl. They showed that by a slight generalization 
of the integrals of motion of the DBCS model, another integrable model is obtained. We shall show that the latter 
can be obtained by taking the quasiclassical limit of the inhomogeneous XXZ vertex model with twisted boundary 
conditions (TIXXZ model), in complete ananology to the derivation of the DBCS model from the TIXXX model. 

Another interesting direction in which our work could be pursued, is to consider boundary conditions with nondi- 
agonal i^T-matrices. These generally lead to models which are not solvable by the ABA. However, their eigenstates 
and eigenvalues can, in many cases, nevertheless be found using Sklyanin's method of separation of variables. 

The paper is intended to be accessible also to readers that are not thoroughly familiar with the details of the algebraic 
Bethe Ansatz; those of its features which are needed here are therefore introduced and reviewed in pedagogical detail. 
The structure of the paper is as follows. In Section || we introduce the DBCS and Gaudin models, recall how their 
integrals of motion are constructed, and give the equations (Richardson's or Gaudin's) that have to be satisfied 



in order to obtain eigenstates and eigenvectors. In Sections |I| and we give a self-contained introduction to 



the ABA method, as applied to the XXX and XXZ vertex models; since both are special cases of the so-called 
6-vertex model, we shall actually begin by discussing the latter in full generality, before specializing later on. In 
Section III we explain how the Yang-Baxter equations satisfied by the R-matrices of local Boltzmann weights lead 
to the exchange relations for the components of the Monodromy matrix T. Furthermore, we derive the fact that the 
transfer matrices for different spectral parameters commute, which is the underlying reason for the integrability of 
the model. In Section IV we exploit the exchange relations of the components of the Monodromy matrix to construct 
the eigenstates and eigenvalues of the model. In Section ^ we explain how the results of Sections III and IV can be 
straightforwardly generalized to the case of twisted boundary conditions using Sklyanin's K- matrix. In Section Vl, 
we take the quasiclassical limit of the TIXXZ model, and show that one thereby recoveres a generalized version of the 
DBCS mode l. We also show that if one specializes these results to the TIXXX model, one recovers the DBCS model. 
Section VII contains some brief conclusions and an outlook for future applications of our results. 



II. THE DBCS AND GAUDIN MODELS 



The DBCS model that is commonly usedat3 to describe superconducting pairing correlations in nanoscale metallic 
ifrains is defined as follows: one consideres a reduced BCS Hamiltonian, 
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^ = ^A.y^O'^ - 9Yc]^c]-C]'-Cj'+ , (1) 

j:0-=± n' 

for electrons in a set of pairs of time-reversed single-particle states |j, ±) with energies Ej, which are scattered pairwise 
from level j' to j, with interaction strength g. Richardson managed to solve this model exactly, for an arbitrary set of 
levels Ej (degenerate levels are allowed, but are to be distinguished by distinct j-labels, i.e. they have Si = Sj for i ^ j). 
Since any level occupied by only a single electron does not participate in and remains "blocked" to the pairscattering 
described by H, the labels of all such single-occupied levels are good quantum numbers. The eigenstates \a) and 
corresponding eigenenergies £a of H thus have the foUwoing general form: 

Here B is the set of singly-occupied, blocked levels, and |-0p) is an eigenstate, with eigenvalue £p and containing 
precisely P pairs of electrons, of a Hamiltonian Hjj which has precisely the same form as the H of Eq. (^, except 
that now the j-sums are restricted to run only over the set U of all unblocked or non-singly-occupied levels. It is now 
convenient to introduce the pseudospin variables 



S. 



i (l - c]+c,+ - c]_c,„) , S+ = c^_c^.+ , = c]+c]_ , (3) 



which satisfy the standard SU{2) relations 

[S+,SJ] = S,, 2 , [S:,Sf] = ±S,,Sf , (4) 
and in terms of which Hir takes the form 

u 

Hu^Y. 2£.(l/2 - - E ^r^; • (5) 
Our choice for the pseudospin variables differs from that used in many other publicationsBi^ by the replacement 

s+ ^ s; , St ^ ~st , (6) 

which preserves the 517(2) relations (||). With our choice, the physical vacuum state |0), containing no pairs, has the 
maximum possible S'z -eigenvalue and hence is a "highest-weight" state. This is convenient for our present purpose, 
namely to establish contact with the ABA, because in the Bethe-Ansatz literature it is standard practice to use 
highest-weight states as reference states [see Eq. ( p6| ) below]. 

Nassj, Richardson showed that the sought-after P-pair eigenstates (unnormalized) and eigenenergies have the general 
forirO 

P (7 

\^p)^\{s-{^lm. with 5-(A.0 = E^r^' 
1=1 t ^' 

p 

Sp = Y.^^l. (8) 

1=1 

Here the P parameters /i/ {I — 1, . . . ,P) are a solution of a set of P coupled algebraic equations, which we shall call 
the "Richardson equations" , 

--T.0 +E =0' 1 = 1,. ...P. (9) 

These are to be solved (numerically, see Appendix B of Ref. |^) subject to the restrictions [ly ^ /i; if V ^ I. A simple 
proof of this result may be found in Appendix B of Ref. |[ its strategy is to verify that [Hu — £p)\^! p) = by simply 
commuting Hjj to the very right past all of tlte operators in (^. l— , 
Moreover, Cambiaggio, Rivas and SaracenoEl showed that the constants of the motion of Hu have the formlla 
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Hi ^ S^+gY^ ^'^^ + UStS-+S-S+) ^^^^ 

The operators Hj, i = 1,2, . . . , N commute with each other as well as with the Hamiltonian (^. In the limit 5 — > oo, 
the operators 

HP-<^'" = hm Hi/5 = y '^^^ + ^(^^^7 + ^^^) (n) 

g^oa ^-^ Ei — En 

j=l ■> 

coincide with the Hamiltonians of the Gaudin chain (see chapter 13 of Ref. The common eigenstates of the Gaudin 
Hamiltonians are given by the same Eqs. (|^), but with the parameters /i/ satisfying the so-called Gaudin equations, 
which are simply the g —> 00 limiting case of Richardsons Eqs. (^). The corresponding eigenvalues of the Hamiltonians 
jjGaudin j^j.g given by 

.Gaudin ^ _ — i — i y (^2) 



III. ALGEBRAIC BETHE ANSATZ FOR THE INHOMOGENEOUS 6- VERTEX MODEL 



A. Definition of Model 



It is well known that the Gaudin model can be obtained by taking the quasiclassical limit of the IXXX model. Q The 
main result of this paper will be to show that a similar construction can be used to obtain the DECS model from the 
TIXXX model, as well as generalized DECS models from the TIXXZ model. To set the scene for these developments, 
the next two sections give a pedagogical review of the AEA as applied to IXXX and IXXZ models. Since both are 
special cases of the so-called 6-vertex model, we begin by discussing the latter in full generality. Eethe Ansatz experts 
may want to skip directly to section ^. 

The 6-vertex model is a classical statistical mechanics model on a two-dimensional regular quadratic lattice, whose 
dynamical variables are arrows living along the horizontal and vertical edges of the lattice, labelled by m = 1, . . . , M 
and i — 1, . . . ,N, respectively. At each vertex, only those 6 configuration of arrows are "allowed", i.e. have nonzero 
local Eoltzman Weights (BW), for which the total flux into the vertex is zero (see Fig.|l|). 

A 1 I A A 1 



^ i i f i I 

a a b b c c 

FIG. 1. The six allowed configurations of arrows at a vertex of the 6-vertex model, with their corresponding Boltzmann 
weigths a, b and c. 



Thus, every allowed configurations has exactly two incoming and two out-coming arrows. Furthermore, we take 
the local EW's to be invariant with respect to the simultaneous reflection of all four arrows. This leaves only three 
independent EW's per vertex, to be denoted by 0^1, and Cmi, where the subscripts give the location (m, i) of the 
vertex (intersection of row m and column i). Since the EW's are allowed to depend on the location of the vertex, we 
are considering an "inhomogeneous" model. As usual, the total statistical weight of any given configuration is defined 
as the product of the EW's of all vertices, and the partition function is defined as the sum of these statistical weights 
over all qallowed configurations. 

It is convenient to associate a two dimensional vector space with each row, say Um for row m, and another with 
each column, say Vi for column i, in such a way that the basis vectors 
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(13a) 



represent right- and left-pointing arrows along row to, or upward and downward arrows along column i, respectively. 
Then the local BWs at vertex (to, i) may be viewed as the matrix elements of a linear operator Rmi that acts as 
follows on the tensor product of the TO-th "horizontal" and z-th "vertical" space, Um 'S) Vi: 

Rmi e|") ® e^:^ = iHmifiK 4"^ ® e^^ , (14a) 

where the usual covention of summation over repeated indices l,k € {1,2} is implied. [As a rule, we shall put bars 
over all repeated indices, and tilde's or nothing over non-repeated upper or lower indices, respectively.] It follows that 
the action of Rmi on the coordinates (wmi)"^ of a general vector Wmi = (wmi)"^ e^™'' ® e^*"* G Um ® Vi takes the form 

(RmiWmi)'^ = (Rmi)'?fc(Wmi)'' . (14b) 

The only nonzero matrix elements of the operator Rmi are 

(Rmi)^!"'"! = (Rmi)^2^2 = CLmi , 

(Rmi)^1^2 = (Rmi)^2"'"l = ^mi , (15) 
(Rmi)^2^1 — (Rmi)^l"'"2 — Cmi ■ 

A convenient matrix representation for Rmi is 



(Rmi)';\ = [ 1 , (16a) 

V (7mi) fe (^mi)\- / ; 



where cxmi, /3mi, /3miJ 7mi a-re operators acting on the two dimensional vertical space Vi'. 

farm 0\ \ /^O c„A A f^"" 0\ 

Even more explicitly, Rmi can be expressed as follows in terms of the unit operator I and the Pauli cr-matrices cr^, 

(T± = {(T^ ±i(Ty)/2, 



(Im ® II) + "'^ (t^ ® ,T?) + (,T+ ® (Xr + ^m ® <T^) , (17) 

where the lower indices of the operators indicate the space (Um or Vi) on which they act. 



B. Monodromy Matrix 

One of the most important objects in the ABA method is the Monodromy matrix T„i- It is defined to be the 
operator 

%n — RmN RmN-1 • • • Rmi , (18a) 

which acts on the space Um ® T^i . . . ® Vat, with each factor Rmi acting nontrivially only on the "horizontal" space Um 
and the "vertical" space Vi. To illustrate this action explicitly, we note that the matrix elements of %n are constructed 
as follows from those of Rmi^ 

fcjv...fel _ /-T} \l kM /T} \'jV-l fejV-1 ^T} \liki l'1aK^ 

i-'™) I, kM..M - y^mN) fc„(,RmN-lj • ■ • (,Rmlj ; k^ ■ U»Oj 

This equation has a simple physical interpretation: each such matrix element of Tm gives the total Boltzmann weight 
of the TO-th row, depicted in Fig.H, for a fixed configuration of external arrows (specified by the indices of Tm on the 
left-hand side of Eq. (|18b|)), obtained by summing over all allowe d con figurations of arrows on internal horizontal 



edges (the k sums, for i = 2, . . . , A^, on the right-hand side of Eq. (18b)). Likewise, using this one row construction 
as building block, the partition function of an AZ-row lattice can be expressed via the matrix elements of a suitable 
product of the M Monodromy matrices, as will be seen below. 



5 



2 1 



m 



W-2 



kj 
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FIG. 2. Construction of the Monodromy matrix: the matrix element {%nf\ ^'^'^'^].^ is equal to the total Boltzmann weight 
of the m-th row, for a fixed configuration of external arrows [specified by the indices of %n\, obtained by summing over all 
allowed configurations of arrows on internal lines (whose indices carry bars here). 



Because of the different roles played by the horizontal space Um (usually reffered as the auxiliary space) and the 
tensor product of remaining vertical spaces Vi ■ ■ ■ Vm (the so called quantum space), it is convenient to arrange the 



matrix elements of 7^ that correspond to the horizontal space C/„ 
2x2 matrix: 



i.e. {Trnfi in the notation of Eq. ( |l8q ), into a 



J. m 



A R 



(19) 



Its entries A, B, C and D are, of co urse, operators acting on the quantum space, which implicitly carry the k- 
indices that are displayed in Eq. ( |l8b| ) (for brevity, we suppressed these above). Each of these four matrix elements 
corresponds to one of four possible kinds of rows in Fig. 2, depending on how the arrows on the first and last (i.e. 
external) horizontal edges are fixed: 



A:{- 



C 



0, D : (. 



(20) 



Note that in this matrix representation, the product of Rmi ma trices on the r.h.s of Eq. (18a) may be viewed 
as conventional multiplication of 2 x 2 matrices [of the form ( 16a )], whose entries are, however, operators on the 
quantum space (and hence carry suppressed /c-indices). 



C. Transfer Matrix 



In order to investigate our model in the case of periodic boundary conditions in the horizontal direction, it is natural 
to consider the operator Tm (called transfer matrix), defined as the trace of %n in the horizontal space Um- 

T„ = Tr„{T„} = J2 C^-^yi = ^™ + ■ (21) 

1=1,2 

Its matrix elements (Tm)'^^^^''^^^ give the total Boltzmann weight of the m-th row depicted in Fig. 2 for a fixed 
configuration of arrows on the vertical edges [specified by the indices of Tm] , obtained by summing over all allowed 
configurations of arrows on horizontal edges, with the boundary condition that the arrows of the first and last 
horizontal edges are equal. 

It follows that the full partition function for a lattice of M rows and N columns can readily be constructed by a 
suitable product of M transfer matrices: for the case of double periodic boundary conditions, it is equal to 

L.M T.M k"'^ k"^^ k^ k^ 

^Af,^ = (.J-Mj rM-l rJ\f-l I, J- M-1 J j-M-2 r Af -2 • ■ • I, J- 1 J rM uM \^^) 

= Tr{TMTM-i...Ti} , (23) 
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i.e. the trace is over the entire quantum space Vi (g) . . . V/y. 

The distinguishing feature of integrable models is that the transfer matrices for different rows commute, TmTm' = 
Tm'Tm- In this case, all transfer matrices have common eigenstates, with eigenvalues Am \ say, so that the partition 
function takes the form 

M 

^A^.^ = E n • (24) 

a m—1 

Thus, the calculation of the partition function reduces to the problem of finding the eigenvalues of the transfer 
matrices. 



D. Yang-Baxter Relations 

It turns out (and will be shown below) that the transfer matrices commute if the local BW's ami, bmi, Cmi are 
parametrized as follows: 

O'mi — 1 ; ^mi = ^(-^miCi) i ^mi — c(Am:^i) ; (25a) 

'^''^^^ ^(A-^ + 2,) ' ^("'^)= ^(A-e + 2,) ' 

where form of the function (j){x) can be either (p^x) — x or 4'{x) = sinhx. The parameter Am (called "spectral 
parameter") is associated with the m-th horizontal line, and (the inhomogeneity parameter) with the i-th vertical 
line. Note that the dependence of the local BW's a^ni, hmi, Cmi on their indices thus enters only via their spectral 
parameters, which is why we could introduce functions 6(A, ^) and c(A, ^) that don't carry the indices (m, i) any 
longer. Note also that the ratio c/b is antisymmetric under an interchange of its arguments: 

c(A,0 ^ _c(e^ 

6(A,0 6(e,A) ' ^ ' 

a property that will be useful later. 

Usually the rational case (j>{x) = x is referred to as the XXX model and the trigonometric case (t>{x) — sinhx 
as the XXZ model, since the Hamiltonians of the XXX and XXZ Hcisenberg magnetic chains can be derived from 
the corresponding (homogeneous) transfer matrices, by taking a logarithmic derivative with respect to the spectral 
parameter at some specific point (se e e.g. chapter 10.14 of Ref. 

For the choice of BW's of Eqs. (p5|), the Rmi-matrices have the following very important property, which ultimately 
leads to the solution of the problem: they satisfies the Yang-Baxter (YB) equation 

■Rriim' (AjTi , A^/ ) Rjni (Am , ) Rm'i ( A^n' ; ) I^m'i ( Am' ; ) -^mi ( Am ; ) -^mm' ( Am ; Ajt^' ) , (27a) 

where the operator products on both sides act on the space Um ® Um' ® Vi, and the arguments in brackets indicate 
explicitly on which parameters the corresponding operators depend. As before, the Rmi operators act on one horizontal 
and one vertical space, Um'SiVi; their nonzero matrix elements are given in Eq. (Esf), with the parameters ami, &mi and 



Cmi as defined in Eq. ( 25a ), with arguments Am and ^i. In contrast, the operator Rmm' acts on two horizontal spaces, 
Um and Um'', apart from this replacement of vertical space Vi by the horizontal space Um', however, the structure 
of Rmm' is exactly the same as that of Rmi^ the nonzero matrix elements of Rmm' are likewise given by Eq. (|l^), 
where now the parameters amm', bmm' and c„im' have arguments A„i and A„i' (i-e. two A's instead of A and ^): 

O^mm' — 1 , bmm' — b(^Xm,^m') , ^mm' — ^(Arn; Am') ■ 

To be explicit, the Yang-Baxter equation implies the following relations between matrix elements of the transfer 
matrices: 

(Rmm')' J Ji (Rmi)' z'^^. (Rm'i)' l'^ k ^ (R-m'i)' Jf^ ^ (R-mi)' j' k (R-mm')' / I' ■ (27b) 

(As usual, we used bars for repeated indices, which are summed over). Graphically, this equation can be represented 
as the equality of the BW's of the configurations depicted in Fig. ||(b), where the left and right diagrams have the 
same configuration of arrows on all external e dges , and sums over all possible configurations of internal indices are 
implied. The verification of the YB equation (E7a) is straightforward (though rather tedious) and reduces to some 
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simple rational or trigonometric identities. (Actual ly, t he solvability of these Yang-Baxter equations dictated the 
choice of parametrization of a, b and c made in Eq. ( pSq ).) As an illustration, the graph in Fig. 3(b) corresponds to 
the following explicit realization of Eq. (|27b|) : 



(28a) 



(R-mm') 2^1 (-E^mi) 2" 2 (-f^m'i) l' 2 ^" (R-mm' ) 1 2 (-Rrni) 2 1 (R-m'i)' 1 2 
= (R-m'i) 2 2 (R-mi) 1 2 (R-mm' ) 2 1 ' 

or, using Eqs. ([T5|), 

which can be verified to hold if ami, ^mi and Cmi have the form specified in Eqs. (Eq). 



(28b) 
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FIG. 3. Graphical depiction of the Yang-Baxter equations: (a): Schematic depiction of the actio n of (R^ 
interchanges the order of the rows m' and m. (b): General graph for the Yang-Baxter equation (27t). (c): Specific graph for 
a particular configuration of external arrows, representing the specific Yang-Baxter equation (p^). Summing over all possible 
configurations of internal arrows consistent with the given choice of external arrows turns out to give two graphs on the left-hand 
side, but only one on the right-hand side. 



Now consider two Monodromy matrices %n and 7^' with identical sets of inhomogeneity parameters ^i, . . . ^at, 

T{Xm) =%n{^m',^l, ■ ■ -^n) , T {Kn') = T-m' {Kn' \ (,1 , ■ ■ ■ ^n) , (29) 

SO that it suffices to display only the functional dependence on the (arbitrary) spectral parameters A,„ or \m' ■ For 
this case, we shall write the components of T(A„i) as A(Am), B{\„i), C{Xm) and D{X,n), and the transfer matrix as 
T(Am), while using Rmm' as shorthand for Rmm'(Am, Am'). A direct consequence of the YB equation is that two 
such Monodromy matrices satisfy the following exchange relations, 

R.mm'^(Am)T(A„i') — T(Am' )T(A,„)Rinm' , (30a) 

or, in terms of matrix elements, 

\l i' ^l, kiv...ki /\ \rj-l' , kN ■■■ki l\ \ _ rj-V , Rn ...ki (\ kff ■■■ki f\ w-p \Tl' (•inU\ 

[ti-min' ) I I' i, kj^...ki^-^"^)-^ I'. fcjv.-.fcA'^™'^ - ^ i'^ kr,..M^^"^'l^ L fc« . . .fci l'^™ A«-mm' j ; • 1>^UDJ 

Fig. 4 is a graphical representation of this relation. To prove it, one successively "pulls" the crossing of the two 
horizontal lines across from the right-most edge of the quantum space to the left-most edge, using the graphical 
representation Fig. 3 of the YB equation ( p7aD . 
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FIG. 4. 
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Illustration of the exchange relations ( pOaD for the Monodromy matrix, (a ) Th e gene ral rel atio n (30a). (b), (c), (d): 
Three specific choices of external arrows, leading to the three Eqs. (p3a), (|33b|) and (33c), respectively. 



Rewriting the Eq. (30a) in the form 

R-mm'^(Am)'^(Am')-^mrn' ~ T {\n')T {\n) , 



(31) 



taking traces over the spaces Um, Um' and using the cyclic property of the trace operation, one immediately concludes 
that the corresponding transfer matrices commute: 



T(A™)T(A™0 - T(A,„OT(A™) 



(32) 



It is the existence of a one-parameter family of commuting transfer matices that makes the exact calculation of their 
eigenvalues and construction of their common eigenstates possible. 



IV. EIGENSTATES AND EIGENVALUES OF THE TRANSFER MATRIX 



Eq.( pOaD represents in a compact form 16 commutation relations among the matrix entries A(Am), B{X„i), C(A„i) 
and D{\m) of the Monodromy matrix T{Xm) (see Eq. (|l9|)). Below we wright down three of them, which are essential 
for solving our eigenvalue eigenstate problem (the full set of relations can be found, e.g. in chapter VII of Ref. p^ ): 



[S(A,„),S(A„,0] =0 



(33a) 
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\ , B{X„,)D{X,^) + ?[^"''^"^i i?(A„)J(A„0 



(33b) 
(33c) 



These three equations correspond to the graphical equations of F ig. 4(b), 4(c) and 4(d) (in that order); for example, 
Fig. 4(d) represents the following specific realization of Eq. (30b) (whose fc-indices we here suppress): 



(-Rrnm' ) 2 1^ 2 ('A^'^)^ 2. ('A^'^' ) ^ (R-mm' ) i 2^ 2 ("A^^i )^ 2 ("A^^iO 



T\ (Am' )T\ (Am ) (R-mm' ) 



2 2 



(34) 
(35) 



Using Eqs. (|15|) and the antisymmetry property (^6|), this readily reduces to Eq. (33c). 

Now we are ready to construct some eigenstates of the model. As reference state we use the following so-called 
vacuum state or highest- weight state. 



(N) 



Vr 



N 



(36) 



which can be visualized as a row of vertical edges, each of which carries an upward-pointing arrow. It is easy to verify 
that 



C(Am)|0) =0, 

A(Am)|0) = |0), 
N 

Dix^m ^Y[b{x,n,^^m. 



(37a) 
(37b) 

(37c) 



To check, e.g., Eq. (37c), recall that the operator D has the graphical representation shown in Fig. 2, with left-pointing 
arrows put on the first and last horizontal edges. The action on |0) implies arranging upward-pointing arrows on all 
the vertical edges above the horizontal line, but then the only allowed arrangement of ar rows on the edges below the 
horizontal line is again a sequence of exclulsively up ward- pointing arrows. By Eq. (18b), this implies that we must 
have N successive 6-type vertices, as specified in Eq. (37c). 

To obtain more general states, one can act on the reference state |0) by an arbitrary number P (with 1 < P < M) 
of operators B: 



\Hi,...,^ip) = B{fj.i)...B{pp)\0). 



(38) 



Note that the parameters /i;, / = 1,...,P, in the arguments of the _B's are arbitrary now, and unrelated to the 
spectral parameters introduced earlier, for which Am is associated with row m. We shall now show that when these 
fii parameters are solutions of a particular system of equations (the famous Bethe-Ansatz equations, see Eq. ( |45| ) 
below), the vector (|8|) becomes an eigenstate of the transfer matrix T(A) (we shall henceforth usually drop the index 
m on Am, since only the functional dependence is important). To this end, let us analyze in detail the action by 
the operator D{X) on the Bethe vector ( p8| ) (subsequently, the action of A{X), which is analogous, will be outlined 
somewhat more briefly). Our strategy is simple: using the exchange relations ( |33c| ), we move the operator D to the 
right past all _B's until it appears next to the vacuum state |0), on which it acts according to (37c). The result can 
be represented as: 



D 



{X)\fii,...,fip) =11 ,./ s n&(A,6)lMi,---,Mp) 

1=1 ^ ' i=l 
p 

+ ^fF\^^l^ ■ ■ • A,/i;+i, ■■■fJ.p) 

1=1 

I'^i 



(39a) 
(39b) 



The first (so-called "wanted") term on the r.h.s. of Eq. (39a) a rises from the case for which one picks up, at each of 
the series of commutations, the first term of the exchange rule (p3q), which in our case is equal to 



10 



b{X,fii 



-B{tii)D{\). 



(40) 



Below we will refer to this term in exchange relation as a "reg ular" term. All other terms (the so called unw anted 
terms), which have contributions from the second term of (33c) , are combined in the second line of Eq. (39a). The 
form of the coefficient fP occuring in this term is given in Eq. (39b) and can be derived as follows: the l.h.s. of the 
Eq. (39a) is symmetric with respect to the permutations of the parameters ni, . . . , fj.p, due to the commutativity 
( |33a| ) of _B's. Since the first term on the r.h.s. of Eq. (39a) is symmetric as well, the second should be too. This means 
that if one succeedes to determine a single coefficient /p for some fixed /, then all other /'s can be straightforwardly 
found using the symmetry. Let us consider the case 1 = 1. It is not difficult to see that the only possibility to obtain 
a t erm that is proportional to |A, /i2, . . . , fip) and docs not contain the operator _B(/ii), is to choose the "wrong" term 
of (|33c|) , 



B{X)D{ni) , 



(41) 



at the very first step when commuting D{\) with i?(/ii), and then everywhere else to choose "regular" ones. Thus, 
for /f we obtain: 



A" 



n 



1 



c{p-i,X) 



N 



(42) 



The above-mentioned symmetry under the permutations of /xj's then immediately implies that in general ff must 
have the form given in Eq. (39b). 

A similar consideration of the action by the operator A on the Bethe vector (5^) gives: 



^(A)|m 



-Aip) = n 



1 



%i,A) 



ti, . . . , 



1=1 



I'^i 



n 



(43a) 
(43b) 



Combining Eqs. (p9) and 



|), we see that the state ( |3^ ) is an eigenstate of the transfer matrix T(A) = A{\) + D{\) 
T(A)|/xi, . . . ,/ip) t(A; ^1,..., /ip) 1^1,..., /ip) (44a) 



with the eigenvalue 



,Mp) =n 



n 



N 



\b{iJ,i,\) fJ^b{X,^ii) fj^ 



(44b) 



provided that + ff = for every I = 1,2, ... P. Using the antisymmetry property of the ratio c/b, Eq. (|2^), it is 
easy to see that this condition is satisfied provided that the P parameters ^i, . . . ,fip satisfy the following system of 
P equations, 



n 



N 



(45) 



I'^i 



which are known as the Bethe equations. Every solution of the system of Bethe equations defines an eigenstate and 
corresponding eigenvalue of the transfermatrix T(A) via (|3q) and (44b). 



V. SKLYANIN'S K-MATRIX 



In this section, we shall generalize, following Sklyanin,lilrlid the formalism described above to the case when the 
boundary conditions in the horizontal direction is not strictly periodic: instead, the I'st and A^' + l'st horizontal bonds 
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are to be identified only up to a "twist", implemented using a (fixed) linear transformation. We will show later on 
that the DBCS model (and also some of its possible generalizations) is some special limiting case of the inhomogeneus 
XXX (XXZljiaodel with such a twisted boundary conditions. Consider a diagonal 2x2 matrix Km, first introduced 



by Sklyanir 



acting on the horizontal space U„ 



(ifm)ll 







{Km)22 



It is easy to check that the following relation holdsQ (illustrated in Fig. |^), 



(46) 



(47) 




(a) (b) 

FIG. 5. (a) Sklyanin's K- matrix and (b) a graphical depiction of Eq. (^). 



where, in accord with our earlier conventions, the subscripts (m, m') specify the horizontal spaces on which the 
operators act non trivially. Let us define a modified Monodromy matrix as 

r„ = r(A„) = K^T{Xm), (48) 

or, m the 2 X 2 block form of Eq. (|l|). 



/ A-m 


Bm \ 


-i 




Dm J 





iKm)22 \Cm D, 



Eq. (47) ensures that the new Monodromy matrix T obeys exactly the same exchange relation (30a) as T, so that 
in particular (|3^) and ( p3| ) remain valid also after the substitution T ^ T, A ^ A, B ^ B, C ^ C, D ^ D. 
Furthermore, the analogues of Eqs. ( |37| ) take the form: 

C(A)|0) =0, 
i(A)|0) =i^n|0), 

N 

D{X)\0)^K22l[b{X,Q\0). (50) 
1=1 

It follows that the Bcthe vector [cf. Eq. (|3|)] 

|/ii,...,/ip)K = ^(Mi)---^(Mp)|0) (51) 

is an eigenstate of T(A), 

f{X)\ni, . . . ,/ip)K = ^(^; /^i: • • ■ ,Mp)Imi, • • ■ ,Mp)k ' (52) 

with eigenvalue 

p p N 

<(A; Ml, . . . , Aip) = i^ii n u + ^22 n uV-^ U ^(^' ' (53) 

provided that the parameters fij satisfy the following Bethe equations, for I = 1, . . . , P: 
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VI. THE "QUASICLASSICAL" LIMIT 

In this section we show how the DBCS pairing model, or a generahsation thereof, can be recovered by taking the 
so-called "quasi-classical" limit (77 — > 0) of the TIXXX or TIXXZ model, respectively. We shall present explicitly 
calculations for the TIXXZ case, i.e. for (j){x) = sinh(a;); to recover the corresponding results for the TIXXX case, one 
simply has to replace all hyperbolic functions by the corresponding rational ones. 



A. Generator for Conserved Operators 



Before taking the limit ry — > 0, it is convenient to write the inhomogeneity parameters as 

^, = 2e, + r], (55) 

where the new parameters e^'s are taken to be independent of rj, and rescale the Rmi operators by a scalar factor, as 
follows: 

R, . ^ 2Rmi, _ sinh(A„, - 2ei + rj) ^ _ 
&(Am,C») + 1 sinh(Am - 2£j) cosh(?7) 

These transformations are convenient because, first, then the leading term in Rmi is simply a direct product of unit 
matrices (see Eq. (57) below); and second, as we will see later, then many equations transform simply under rj —rj 
(being either symmetric or antisymmetric), which considerably simplifies all expansions in powers of r/. When written 
in terms of a direct product of 2 x 2 Pauh matrices as in Eq. ([l7|), the rescaled Rmi of (Eq. (^) takes the form 

■r^ T T tanh?7 ^ ^ 2sinh?7 / , _ is ,^„, 

Rmi = Im ® li + — ^TT ^<^m ® '^f + ^l^T ^ G^m ® '^i + '^m ® '^i^ • 57 

tanh(A„i - 2ei) smh(Am-2£i)' 
Now choose the following form for the K-matrix: 

K = I+^(t\ (58) 

g 

and expand the transfer matrix in powers of the parameter rj, using ( ^7| ) and (|58|). This readily yields 

f (A„,) = Tr„,{KmRmN ... Rmi} = 21 + -V^P{Kn) + 0(7?'), (59) 
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where 



1 ^ 



2 ^ tanh(A - 2e,) 

. ^ / i j i j i j \ 

I 2tanh(A - 2ei) tanh(A - 2ej) ^ sinh(A - 2e^) sinh(A - 2ej) j ' 

i<j 

and to the commutativity (^2|) of transfer matrices for different spectral parameters guarantees that 

[P(A),P(A')]-0. (61) 

P(A) can be viewed as the generating operator for all possible conserved operators of the model. 

A convenient way of obtaining a complete set of commuting conserved operators, is to take the residues of P(A) at 
the points A — 2ei for i = 1 . . . iV, 

Hi = Rcs[P(A); A ^ 2e,\ = I ^P{X) , (62) 

where Ci is a small contour in the complex A-plane, encircling the point 2si. Explicitly evaluating the residues for the 
present model, one obtains 



13 



erf ^/ a? a? (T^crJ- +(Tr(Tf\ 

~ T 1^2tanh(2ei - 2ej) ^ sinh(2e,; - 2ej) j ' ''^^^ 

Eq. ( |6ll ) immediately implies that all of these operators commute: [Hi,Hj] = 0. Furthermore, it is not difficult to 
show that the set of all Hi is complete, in the sense that P(A) can be expressed purely in terms of these operators. 
Indeed, P(A) is a rational (matrix- valued) function of the variable z = exp(2A), which is regular at z — > oo and has 
simple poles at z — > exp(4£i); it is thus completely determined in terms of the corresponding residues, which are equal 
to 2exp(4ei)Hi, so that we have 

P(A) = P(oo)+f:J^. (64) 
The term P(oo) itself also can be expressed via Hi: 

P(oo)=g(f]Hi) +f^Hi-^, (65) 



where we have used the fact that 

N ^ N 

2 



N N 

EHi = ;jE-?- (66) 



The commuting operators Hj are in fact just the so-called generalised Gaudin Hamiltonians. Moroever, in the limit 
g —>■ cxD, in which K = I so that one recovers periodic boundary conditions, the Hi/g reduce to the standard Gaudin 
Hamiltonians Hp^"'^'" of Eq. (|ll|). 



B. Eigenvectors 



To obtain the quasiclassical limit of the Bethe eigenvectors . . . ,/ip)K °f ^'i- (|5l|), we have to investigate the 
77 — > limit of the operators B which enter in its definition of (pT|). Now, recall that the operator B,n = -B(A„i) is the 
(1,2) component (in auxiliary space Um) of the Monodromy matrix 7^ = KmR-mN • • ■ R-mi, which has the following 
expansion [using Eqs. (|57|) and (|5S 



f(A„0 = I+-'T 

g 

N 



(67) 



2?? 



tanh(Am - 2ei) 



sinh(An 



2e^ 



(cr+CT; +(T^cr^ 



In this equation, the only terms having non-zero (1,2) components and hence contributing to B[\m) are those 
proportional to cr+j. Thus we have 



i?(A) = 27?5-(A) + 0(772), 



where 



N 



^ sinh(A — 2ei) 

We see that the quasiclassical (unnormalized) Bethe vector of Eq. (js^) takes the form 

\Hl...^^.p)^ = S~{^ii)...S^{lip)\{)), 
where, as before, the reference state |0) is defined by (^). 



(68) 



(69) 



(70) 
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C. Eigenvalues 



The eigenvalues hi of the conserved operators Hi can be found from the eigenvalue p{X) of their generator P(A). 
Since (4/g)P(A) is the order-77^ coefficient of the transfer matrix T(A) [Eq. (^)], its eigenvalue (4/g)p(A) is given by 
the order-?7^ coefficient of the corresponding eigenvalue t{X), which can be found by multiplying Eq. (p3) by the factor 



N 

n 



sinh(A — 2ei + ij) 
L sinh(A — 2ei) cosh(?7) 



(71) 



[cf. (]5^)] and setting = 2e^ + ?; [cf. (|55|)] 

N 



f(A) 



1 



?/ \ y-r sinh(A - 2ei + r/) -p-|- sinh(A - fJ^i — 2ri) 
g J sinh(A — 2ei) cosh(?7) j-^ sinh(A — /i;) 



(72) 



N 

n 



sinh(A — 2ei — rj) 



n 



sinh(A — /i; + 277) 
9 J sinh(A — 2ei) cosh(?7) J-^ sinh(A — fii) 

Expanding this expression in 77, the coefficient of 77^, multiplied by g/4, is found to be: 

1 ^ 1 



1 ^ 



2 ^—^ tanh(A — 2ei) ^—^ tanh(A — /i/) 
p 



2.9 



1 9 

2 ^-—^ tanh(A — 2ei) tanh(A — 2ei') ^ tanh(A — jii) tanh(A — /i;/) 



KV 
N P 

■EE 



i I 



tanh(A — 2ei) tanh(A — /i/) 



(73) 



The eigenvalues of the generalized Gaudin Hamiltonians Hi, say hi, can be obtained by taking the residues at the 
points A = 2ei: 



1 ^ 
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2 tanh(2e., - /i;) 2 ^ tanh(2£.j - 2ei' 



(74) 



Finally, since all the Hi commute, we can immediately write down the eigenvalue of any function of these operators. 
In particular, the general Hamiltonian 



H = P(Hi,...,Hn), 



(75) 



where V is some arbitrary polinomial of its arguments, has eigenvalues V{hi, . . . /iat). For example, the general class 
of models recently discussed by Amico, Di Lorenzo and Osterlohu in the context of superconductivity in small grains, 
is obtained by considering certain second order polynomials (i.e. quadratic combinations of Hi's). 



D. Bethe equations 



The quasiclassical Bethe state . . . ,/ip)K is an eigenstate of the generator P(A), and consequently also of each 
of the generalized Gaudin Hamiltonians Hi, only if the parameters /i/ satisfy the limit 77 of the Bethe equations 
(^). The latter are of course not affected by the rescaling transformation (^6|), and take the following form upon 
inserting ^i = 22^ + 77 of ( ^ ) : 

(l+'l\r\ sinli(Mi' - Mi + 27?) ^ f^ _v\-rT sinh(^; - 2ei - 77) 
V g J j;^ sinh(^z' - fii - 2rj) \ g J fj^ sinh(^/ - 2ei + 77) 
I'^i 



15 



In the "quasiclassical" limit ?/ ^ we obtain the following set of equations, for I = 1, .... P, which may be viewed as 
generalized Gaudin equations: 



N 



g ^ tanh(2ei 



E 

i'=i 



tanh(/i;/ - /i;) 



= 0. 



(76) 



We would like to emphasize that these are the on-shell Bethe equations of quasiclassical limit of the TlXXZ model. 
In contrast, in Rcfs. ^ and |9| who did not consider twisted boundary conditions as we do here, these equations are 
off-shell Bethc Ansatz equations. 

Of course, Eqs. (|76|) can be derived, if desired, without reference to the ABA, by pursuing the following strategy 
(described in detail in Appendix B of Ref. 3|): in order to show that the state . . . Mp)k -'^l- ( [^ ^-i^ eigenstate 
of any Hi, one would commute Hi past all the operators S~{fii) in Eq. ( |7C| ) [whose form ( |69| ) is reminiscent of the 
operators defined Ref. || if we identify o-r with 5^]; this would generate "unwanted" terms that only vanish if 
Eqs. (|76|) are satisfied. 



E. Specialization to Richardson's equations 

It is straightforward to recover the DBCS model and Richardson's solution thereof, as summarized in Section ||, 
by considering the case (t>{x) = x appropriate for the XXX model (instead of the XXZ case = sinhx), and replacing 
everywhere 



tanh; 



sinh: 



(77) 



First, we note that the generalized Gaudin equations ( |76| ) then reduce to Richardson's equations (g). Furthermore, 
the generalized Gaudin Hamiltonians Hi of ( pS] ) reduce to the form given in Eq. ( |l0| ) for the conserved operators 
of the DBCS model. This fact was noted by Sklyanin himself in a side remark in Ref. |ll], and first derived by him 
already in 1989 in Ref. |l^. However, he was at the time unaware of the fact that the resulting H^ were useful in 
the context of the DBCS model, and in particular, that they be used to construct the Hamiltonian Hu of (H) of the 
DBCS model. It it is straightforward to check that this can be done through the following construction: 



N N 

ifE7(Hi)=^[(g-2e,)Hi + (£,-35/4)] + 9{Y.^^ 



(78) 



To calculate its eigenvalues £p = Hij{hi) explicitly, the following identities [derived by repeated use of (74) and (|7^)] 
are useful: 



N 



N 



N 

N 



P N 

EE 

1=1 i=l 
N 



2ei - 



N 



9^1 



N 



7^' 



E 

1=1 

p 

EE 



- ( 

i+E- 

l'=l 

_2ge, 
2f 

1=1 1=1 
p 



2.9 



N 



P, 



N P 



2£i - ^/ 

^ + gN[N - l)/4 - A^/ + qPN - gP{P - 1) 



(79a) 



(79b) 



(79c) 



i=l 



1=1 



where the last term of (79c) was obtained from the last term of (79b) by rewriting the latter as follows 



N P 



i=l 1=1 1=1 i' = i 



2.9 



(79d) 



Using Eqs. (79a) and (79c), it is straightforward to check that £p reduces to the simple form £p ~ J2i=i M' of Eq. 
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VII. CONCLUSIONS 



The TIXXZ results of the previous section for the conserved operators Hi, their eigenvalues hi and eigenvectors 
1^1, ... , /^p)k) ^'iid corresponding consistency condition ([76|), have been found independently before by Amico, Di 
Lorenzo and OsterlolJj. They managed to construct the Hi apparently by inspection, without presenting a systematic 
approach for their derivation, and in their approach the consistency condition (|7^) appears as a set of off-shell Bethe 
Ansatz equations. In our work, we presented a systematic derivation of these results from a vertex model with twisted 
boundary conditions, and the consistency condition ( [76| ) corresponds directly to the on-shell Bethe Ansatz equations 
of this model. Thus, we hope to have have shed some additional light on the reasons why the DECS model and its 
generalizations are integrable and Bethe- Ansatz solvable, and on the underlying algebraic structure of the solutions. 
We hope that our works shows the way towards further progress in applyins-the powerful formalism of the ABA to 
the DBCS and related models, e.g. for the calculation of correlation functionalj such as (SfSj) or {Sj^S^l. which are 
of importance for understanding the nature of pairing correlations in nanoscale superconducting grains .013 

Acknowledgements: — We acknowledge very fruitful discussions with and stimulating encouragement from R. Flume 
and V. Rittenberg, and helpful comments from L. Amico. RP has been supported in part by the INTAS grant 99-1459 
and the A. v. Humboldt foundation of the Federal Republic of Germany. 



D.C. Ralph, C.T. Black and M. Tinkham, Phys. Rev. Lett. 76, 
' J. von Delft, Annalen der Phys. (Leipzig) 10, 219 (2001). 



(1996); 78, 4087 (1997). 



J. von Delft and D. C. Ralph, Phys. Rep. 345, 61 (2001). 
* R.W. Richardson, J. Math. Phys. 18, 1802 (1977), contains a complete set of references on the subject, starting from R.W. 

Richardson, Phys. Lett. 3, 277 (1963). 
^ M. C. Cambiaggio, A. M. F. Riv as. and M. Saraceno , Nucl. Phys. A 624, 157 (1997). 



cond-mat/0010349 



L. Amico, G. Falci and R. Fazio, 
^ M. Gaudin, La fonction d'onde de Bethe; MASSON, 1983. 
® H.M. Babujian, in Proc. XXIV Int. Symp. Ahrenshoop, Gosen 1990; H. M. Babujian, J. Phys. A 26, 6981 (1993), J. Physics 

A 27, 7753 (1994); H. M. Babujian and R. Flume, Mod. Phys. Lett. A9, 202 9 (1994). 

^ L. Amico, A. Di Lorenzo, and A. Osterloh, Phys. Rev. Lett. 86, 5759 (2001); |cond-mat/0105537 . 
E. K. Sklyanin, in "Quantum Group and Quantum Integr able Systems" (Nankai Lectures in Mathematical Physics), Ed. 

M.-L. Ge (World Scientific, 1992) p. 63-97 lhep-th/92imi[l, 

" E.K. Sklyanin, Progr. Theor. Phys. Suppl. 118, 35 (1995) |iolv-int/9504001 1. 
E. K. Sklyanin, J. Sov. Math. 47, 2473 (1989). 

After the bulk of this work had been c ompleted, a preprin t appeared in which some such correlation functions indeed were 
calculated: L. Amico and A. Osterloh, cond-mat/0105141. 

It was recently pointed out by Sierra that the form of Richardson's wavefunctions (Q) can be related to that of certain 
correlators of the su{2)k Wess-Zumino-Witten conformal field theory, in the singular limit when the level k goes to —2, see 
G. Sierra, Nucl. Phys. B572, 517 (2000). 
^® The form for Hi given in Ref. ^ is related to our Eq. (|l^) by the transformation 
R. J. Baxter, Exactly Solved Models in Statistical Mechanics: Academic Press, 1982. 

V. E. Korepin, N. M. Bogoliubov, A. G. Izergin, Quantum Inverse Scattering Method and Correlation Functions; Cambrige 
University Press 1993. 

In fact, it is possible to satisfy Eq. (^), and hence maintain integrability, using more general choices of K: for the XXX 
model, Eq. (^ holds for an arbitrary 2x2 K-matrix. For the XXZ model, mm holds for any diagonal K-matrix, or any 
purely-off-diagonal K-matrix (i.e. vanishing diagonal elements). For both the XXX and XXZ models, though, only the case 
of a diagonal K-matrix is solvable by conventional Bethe- Ansatz methods. The off-diagonal cases are solvable, however, using 
Sklyanin's method of separation of variables. 



17 



